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Abstract 

We compute the imaginary part of the heavy quark contribution to the photon polarization 
tensor, i.e. the quarkonium spectral function in the vector channel, at next-to-leading order 
in thermal QCD. Matching our result, which is valid sufficiently far away from the two-quark 
threshold, with a previously determined resummed expression, which is valid close to the 
threshold, we obtain a phenomenological estimate for the spectral function valid for all non- 
zero energies. In particular, the new expression allows to fix the overall normalization of the 
previous resummed one. Our result may be helpful for lattice reconstructions of the spectral 
function (near the continuum limit), which necessitate its high energy behaviour as input, 
and can in principle also be compared with the dilepton production rate measured in heavy 
ion collision experiments. In an appendix analogous results are given for the scalar channel. 
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1. Introduction 



Heavy fermion vacuum polarization, i.e. the contribution of a massive fermion species to 
the (imaginary part of the) photon polarization tensor, or to the spectral function of the 
electromagnetic current, is one of the classic observables of relativistic quantum field theory: 
the result has been known up to 2-loop, or next-to-leading, or 0{a s ) level already since the 
1950s [T]0 Nevertheless, significant new insights were still obtained in the 1970s [2] and 
even in the 1990s [3]. By now a lot of information is also available concerning corrections 
of 0(ot1) and 0(a^) (for recent work and references, see ref. [4]). The physics motivation 
for the continued interest is related, for example, to determining the heavy quark production 
cross section, a(e~e + — ► cc), often expressed through the i?-ratio, as well as to computing 
the heavy quarkonium decay width. 

In the present paper, we consider essentially the same observable as in the classic works, 
but in a situation where the heavy quarks live at a finite temperature, T, rather than in 
the vacuum. We refer to this observable as the "heavy quark medium polarization" . Again 
the result has direct physical significance, in that it determines the heavy quark contribution 
to the production rate of lepton-antilepton pairs from the thermal plasma [5] (cf. eq. (|2,2p 
below). There has been considerable phenomenological interest particularly in what a finite 
temperature does to the resonance peaks in the spectral function, given that this might yield 
a gauge for the formation of a deconfined partonic medium [6]. Some recent work on the 
resonance region within the weak-coupling expansion, taking steps towards a systematic use 
of effective field theory techniques to resum appropriate classes of higher loop orders, can 
be found in refs. [7]-|llj (see ref. [12j for a review and ref. [13] for an alternative approach 
with similar results), and recent reviews on some of the phenomenological approaches on the 
market can be found in refs. |14t [T5] (possible pitfalls of ad hoc potential models at finite 
temperatures have been reviewed in ref. [16J, and underlined from a different perspective 
in ref. [13]). Analogous spectral functions can also be determined for theories with gravity 
duals J?]. 

Unfortunately, it appears that ultimately weak-coupling (and related) techniques will be 
insufficient for determining quantitatively the shape of the spectral function around the res- 
onance region. The reason is that field theory at finite temperatures suffers from infrared 
problems, implying that the weak-coupling series goes in powers of (a. s /ii) l l 2 rather than 
a s /ir, often with large (sometimes non-perturbative [18]) coefficients; see, e.g., ref. [19] and 
references therein. Therefore, particularly for the case of charmonium where even at zero tem- 
perature weak-coupling computations can hardly be trusted, it appears that non-perturbative 
techniques are a must. Even though the situation should be somewhat better under control 

1 The older computations were formulated within QED, but at this order the results carry over directly to 
QCD, whose notation we adopt. 
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for bottomonium, a crosscheck by lattice methods would still be more than welcome. 

On the point of lattice techniques, however, we are faced with a rather fundamental prob- 
lem. Lattice techniques are applicable in Euclidean spacetime, while the spectral function is 
an inherently Minkowskian object. In principle the spectral function can be obtained through 
a certain analytic continuation of the Euclidean correlator; however, as a mathematical op- 
eration, at finite temperatures such an analytic continuation is unique only if the asymptotic 
behaviour for large Minkowskian arguments is known (see, e.g., ref. [20]). In a practical 
setting, many further problems arise because lattice data is not analytic in nature; yet the 
need to input outside information ( "priors" ) to the analysis certainly remains a central issue 
(see refs. [21] for recent lattice results, and ref. [22] for an overview). 

It is at this point that weak-coupling techniques may again become helpful. The goal 
would now be not so much to determine the spectral function around the resonance region, 
but to determine it at very high energies, which information should be relatively reliable, 
thanks to asymptotic freedom. Indeed, the free thermal spectral function has been studied 
in great detail previously, even at a finite lattice spacing [23]. Nevertheless, loop corrections 
are expected to remain substantial even up to energy scales of several tens of GeV, so it is 
important to account for them, and this is the basic goal of the present work. As a more 
refined goal, we wish to demonstrate that thermal corrections are small far away from the 
threshold; thereby knowledge of the asymptotic behaviour could be taken to higher orders, 
by employing well-known numerically-implemented results from zero temperature [24] (this 
assumes, of course, that a continuum extrapolation can be carried out on the lattice). 

Apart from this lattice-related goal, we also wish to pursue the complementary goal of 
treating the bottomonium spectral function without any exposure to the often hard-to-control 
systematic uncertainties of lattice simulations. This can be achieved by constructing an 
interpolation between the asymptotic result determined in the present paper, and the near- 
threshold behaviour estimated within a resummed framework in ref. [8]. 

The plan of the paper is the following. In sec. [2] we define the observable to be computed. 
The general strategy of the computation is discussed in sec. and the main results are 
summarized in sec. HI A phenomenological reconstruction of the spectral function in the 
whole energy range is carried out in sec. [5j while sec. [6] lists our conclusions. In appendix 
A, we display in some detail the intermediate steps entering the determination of one of 
the "master" sum-integrals appearing in the computation; in appendix B, we list the final 
results for all the master sum- integrals; and in appendix C we provide results for the spectral 
function in the scalar channel, discussing briefly also the ambiguities that hamper this case. 
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2. Basic definitions 



The heavy quark contribution to the spectral function of the electromagnetic current can be 
defined as 

p v (u) = j^dte^ jd 3 - 2 ^ ^[^(t,x),^(0,O)]\ , (2.1) 

where = %[) ^ ip; ip is the heavy quark field operator in the Heisenberg picture; (...) = 
iJ~ 1 Tr [(...)e~^ ] is the thermal expectation value; (3 = 1/T is the inverse temperature; and 

we assume the metric convention (H ). This spectral function determines the production 

rate of muon-antimuon pairs from the system [5] , 

djV-„+ -2e 4 Z 2 / 2ml\ ( 4m 2 A ^ 

where Z is the heavy quark electric charge in units of e, and n-Q is the Bose-Einstein distri- 
bution function. In defining eq. (|2.1|) and the argument of py in eq. (|2.2|) . we have assumed 
that the muon-antimuon pair is at rest with respect to the thermal medium, i.e. q = (u, 0)|§ 
The pole mass of the heavy quark (charm, bottom) is denoted by M. 

The parametric temperature range we concentrate on in this paper is the one where the 
"quarkonium" resonance peak disappears from the spectral function py [8]: 

g 2 M <T <gM . (2.3) 

This implies that in any case T <C M, so that exponentially small corrections, ~ exp(— /3M), 
can well be omitted. The thermal effects come thereby exclusively from the gluonic sector, 
where no exponential suppression takes place. 

In order to compute the spectral function py of eq. (|2,ip . we start by determining the 
corresponding Euclidean correlator, 

C E (oj n ) = J dre^J d 3 - 2e x (^(r !X )jf M (0,O)) , (2.4) 

for which a regular path-integral expression can be given (i.e., hats can be removed from the 
definition). Here ui n = 2irnT, n E Z, denotes bosonic Matsubara frequencies. The spectral 
function is then given by the discontinuity (see, e.g., refs. [25} 126]) 

p v (u) = Disc \c E (-iuj)\ = — \c E (-i[u + i0 + }) - C E (-i[u - i0 + ])] . (2.5) 
L J 2i L 

In the following we denote Euclidean four-momenta with capital letters, in particular Q = 
(u n ,0). Moreover, f K = T^ kn M 2e / d d k/(2vr) d stands for a sum-integral over bosonic Mat- 
subara four-momenta, while $ {P} signifies a sum-integral over fermionic ones. The space-time 
dimensionality is denoted by D = 4 — 2e, and the space dimensionality by d = 3 — 2e. 



2 For a non-zero total spatial momentum q, with < |q| <C M, the main modification of our results would 
be a shift of the two-particle threshold from u) « 2Af to uj ~ 2M + q 2 /4M. 
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3. Details of the computation 



3.1. Propagators 

At a finite temperature T it is not clear, a priori, whether the result of the computation will 
be infrared finite, given that (after analytic continuation) the gluon propagator contains the 
Bose-enhanced factor n B (k°) ~ T/k°, for \k°\ <C T. For this reason, we carry out the analysis 
by using the Hard Thermal Loop resummed [27\ [28] form of the gluon propagator, which 
takes into account Debye screening, and thereby shields (part of) the infrared divergences. 
Introducing (see, e.g., refs. [25} [26]) 

PZ (K) = P£(K) = P£(K) ee , P?(K) ee Sij - ^ , (3.1) 



P^{K) ee (3-2) 



where K = (/c n ,k), k n = 2imT, the Euclidean gluon propagator can be written as 
(A a Jx)A b v (y)) = 5 ab i e iK <*-ti 



A 

T 



PU K ) , P ^( K ) tK^Ky 
K 2 + U T {K) K 2 + U E (K) (K 2 ) 2 



(3.3) 



where £ is a gauge parameter. The projector P is transverse both with respect to K and to 
the four-velocity of the heat bath and, in the static limit, describes colour-magnetic modes; 
the projector P E is transverse only with respect to K and, in the static limit, describes 
colour-electric modes. The self-energies T1t,T1e are well-known \27\ 128] functions of the form 
m D /(^n/|k|), where ttid = (N c /3 + N[/Q)gT is the Debye mass parameter; we will not need 
their explicit expressions in the following, apart from knowing that / is an even function of 
its argument and regular on the real axis. The fermion propagator has the free form, 

W*tf(y)> = i e^~^ ~%+™ 2 ' (3 ' 4) 
{P} B 

where Mb is the bare heavy quark mass. 
3.2. Contractions 

The first step of the computation is to carry out the Wick contractions and the Dirac 
traces. At 1-loop level, omitting Q-independent terms which are killed by the discontinu- 
ity in eq. (|2.5p . we get 

O-- H - W «^ ffi^ff . ,3,) 
Here Ca = N c , and 

A(P) = P 2 + M 2 . (3.6) 
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At next-to-leading order (NLO), we have to evaluate the counterterm graph as well as 
genuine 2-loop graphs. The counterterm graph can be deduced from the 1-loop expression in 
eq. (|3.5p . by re- interpreting the mass parameter as the bare one, JW|, and then expanding it 
in terms of the pole mass: 



M 



,.2 6g CpM fl h 2 4 . 



(3.7) 



where C F = (N 2 - l)/2iV c , and is the scale parameter of the MS scheme. This yields 



[Q - indep.] + 



24g 2 C A C F M 2 (\ fl 2 4 



{P} 



(4tt) 2 

(D - 2)Q 2 - AM 2 2 
A 2 (P)A(P-Q) + A(P)A(P-Q) 



(3.8) 



For the genuine 2-loop graphs, we make use of the identities 



K fM pT{K) = Q tl pT(K) = 0, P^(K)=D-2, P^P^K) = p 2 - (p • k) 2 , (3.9) 

where k = k/|k|, and the second equality follows from the fact that Q is aligned with the 
heat bath. We then complete squares in the numerator, and note that 

QK 



K{P} 

as can be shown with the shifts P 



[K 2 + U(K)]A(P)A(P - Q)A(P - K)A(P - Q - K) ° ' 
P + Q, K — > —K . Thereby we arrive at 



(3.10) 



<^^> + ={q- inde p-] + ^ 2 c A c F i { 

T K{P} 1 

1 1 



(3.11) 



[p 2 - ( P • k) 2 ] x 



K 2 + U T K 2 + U Er 

4[(2-D)Q 2 + 4M 2 ] 



2[(2 - D)Q 2 + 4M 2 ] + 4fT 2 



A 2 (P)A(P - Q)A(P - K) A(P)A(P - Q)A(P - K)A(P - Q - K) 



D 



K 2 + U T 



+ 



+ 



(2 - D)Q 2 + 4M 2 



A(P)A(P-Q) A 2 (P)A(P-Q) 

2 -2(D - 2)Q ■ K + 4K 2 



+ 



A(P)A(P - Q - K) A(P)A(P - Q)A(P - K) 

[(2 - D)Q 2 + 4M 2 ]K 2 [(6 - D)Q 2 /2 + 2M 2 ]K 2 + K 4 



A 2 (P)A(P - Q)A(P - K) A(P)A(P - Q)A(P - K)A(P - Q - K) 



+ 



K 2 + U E 



4[(2 - D)Q 2 + 4M 2 



+ 



4[(2 - D)Q 2 + 4M 2 ]M 2 



A(P) A(P - Q) A(P - K) A 2 (P)A(P - Q)A(P - K) 

(2 - D)Q A + (8 - 2D)Q 2 M 2 + 8M 4 + [(2 - D)Q 2 + 4M 2 ]K 2 
+ A(P)A(P - Q) A(P - K)A{P -Q-K) 
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Note that any dependence on the gauge parameter £ has disappeared; thus P^ v could have 
been replaced with 5^ — P^ v . 

3.3. Outline of the subsequent steps 

Given eq. (|3.1ip . we need to carry out the Matsubara sums and the spatial momentum inte- 
grals. More concretely, the steps (specified in explicit detail for one example in appendix A) 
are as follows: 

• Writing the gluon propagator in a spectral representation, the Matsubara sums T J2k n 
and T^2^ pn y can be carried out exactly in all cases. 

• The result after these steps contains many appearances of the Fermi distributions, 
n F (E) = l/[exp(f3E) + 1], where the energy E is that of a heavy quark, E > M. All 
such terms are suppressed by at least e~ M l T <c 1, and can be omitted. 

• The remaining temperature dependence appears as Bose distributions with the gluon 
energy, n B (k°). Here the issue is the opposite: in the small energy range, \k°\ <C T, 
there is an enhancement factor T/k°, which could lead to infrared divergences. This is 
an important point, so we devote a separate subsection to it (sec. I3.4p . The upshot is 
that there are no infrared divergences at the present order. 

• Having verified the absence of infrared divergences, we can forget about the HTL re- 
summation in the gluon propagators, i.e. set Hr = = in eq. (|3.1ip . and insert 
the free spectral function for the gluons. Thereby the integral over the gluon energy 
k° is trivially carried out. (In practice, we first insert the free gluon spectral function, 
integrate over k°, and verify the absence of infrared divergences a posteriori for each 
independent ("master") sum-integral separately.) 

• The remaining spatial integrals, over k and p, are effectively three-dimensional (over 
the absolute values of k, p and over the angle between k and p). Some of them are 
ultraviolet divergent, and require regularization. The integrals come in two forms, which 
we call "phase space" and "factorized" . We are able to carry out two of the integrations 
in all cases; for the zero-temperature parts entering the final result, all three integrations 
are doable pQ-[3], while for the finite-temperature parts an exponentially convergent 
integral over k = |k| remains to be carried out numerically. 

The results obtained after these steps are listed for all the master sum-integrals appearing in 
eq. (|3.1ip in appendix B. 
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3.4. Absence of infrared divergences 

Inserting the free gluon spectral function, which sets kP = k, into any of the master sum- 
integrals, there remains an integral over the gluon momentum k to be carried out. In principle 
this integral could be infrared divergent. This turns out indeed to be the case for the "phase 
space" and "factorized" parts (for definitions, see appendix \A.2\i of the integrals separately; 
in fact, the integrals denoted by S® and Sq (cf. eqs. ()B.20p . (|B.32p ). have logarithmically 
divergent infrared parts even at zero temperature, which were an issue in the 1970s [2]. 
However, the infrared divergences were found to cancel in the sum of the phase space and 
factorized parts. In our case, the logarithmic divergences turn into linear ones, due to the 
additional factor n B (k) T/k; nevertheless, when the phase space and factorized parts are 
added together, we find that both powerlike and logarithmic divergences cancel, and the 
integrals become finite, for each master sum-integral separately. This can clearly be seen 
m eqs. (IB~24) and (jR35j) for Sg and Sg, respectively. The same is true for the integrals 
denoted by S$,S$,S% (eqs. (jR25]) . (frTMj) . (fTOKIl ). appearing in the first term of eq. (13TTT1) 
and disappearing if the HTL self-energies are set to zero from the outset. Therefore, we 
conclude that there are no infrared problems in our observable at the next-to-leading order. 
(It is to be expected, though, that there are some at higher orders.) 



4. Final result 



Given the considerations in sec. 13.41 showing the absence of infrared divergences, we are free 
to set Ht = IIe = in eq. (|3.1ip . Noting furthermore that the factorized gluon tadpole reads 



J K 



1 T 2 

5, = ia+°M. 



(4.1) 



and employing the notation of appendix B for the sum-integrals Sf (w), the full result can be 
written as 

M»)L = -"Cm^ 2 + 2M 2 )SiM + Sg 2 c A c F 

T " + + 1)} [- SlM + ^ + 2M2 " e ^ (w) 

+2S 3 (uj) - 4(w 2 + 2M 2 - euj 2 )Sl{uj) - 4(1 - 2e)S\{u) + 4(1 - e)Sj(oj) 
+2{uj 2 + 2M 2 - ecv 2 ) \2M 2 Sl{uj) - (1 - e)Sf - (w 4 - 4M 4 )5 6 °(w) 



+ 



(2 - e)uj 2 + 2eM 2 S|(w) - (1 - e)S 4 (w) } + 0(e) . 



(4.2) 



We have set here e — ► whenever the master sum-integral that it multiplies is finite. 
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Now, the explicit thermal correction on the second line of eq. (|4.2p has a simple physical 
meaning: it corresponds to an expansion of the leading-order result through a thermal mass 
shift [29] 

g 2 T 2 C F 



M 2 + 



6 



(4.3) 



i.e. SM = g 2 T 2 C F /12M. Note that this term multiplies the function S 2 (u) = 0(uj - 
2M)/[16ituj(uj 2 - 4M 2 ) 1/2 ](1 + 0(e,e~P M )) (cf. eq. jEfiJ)), which diverges at the threshold, 
while the sum of all the other terms turns out to remain finite. Thereby the thermal cor- 
rection would completely dominate the result close enough to the threshold, were it not to 
be resummed into a mass correction a la eq. (|4.3p . On the other hand, once it has been re- 
summed, this term is in general small: in the range that we are interested in, g 2 M < T < gM, 
it corresponds parametrically to a higher order contribution. Therefore, for simplicity, we 
drop this term in the following (of course, if desired, it is trivial to include it as an overall 
mass shift), and reinterpret the result as 



-AC a (uj 2 + 2M 2 )S 1 (uj) + 8g 2 C A C F 



6M Z 



(47T) 



1 , V 4 
- + In -£=r + - 



M 2 



■Si(u) + (uj 2 + 2M 2 - eu 2 )S 2 (u) 

+2S 3 (uj) - A(uj 2 + 2M 2 - eu?)Sl{uj) - 4(1 - 2e)S\{uj) + 4(1 - e)Sj(uj) 
+2(uj 2 + 2M 2 - ecv 2 ) \2M 2 Si(u) - (1 - e)5f - (uj 4 - 4M 4 )5^(w) 



+ 



[2 - e)u 2 + 2eM 2 



S 2 6 (u)-(l-e)Sl(u) \+0{e). 



(4.4) 



Nevertheless, it is perhaps appropriate to stress that only the part of the thermal correction 
multiplying the function Sz(<jj) can be unambiguously resummed on the grounds that the 
result would otherwise diverge at the threshold, while the term ~ T 2 Si(u) could in principle 
be kept explicit, and would then have an 0(1) effect on the thermal part of the result. 

Inserting the explicit expressions for the functions SHuS) from appendix B into eq. (|4.4j) ■ 
the final result for the vacuum part becomes 



Pv{u)[ 



-9{uj - 2M) 



C a (uj 2 - 4M 2 p (uj 2 + 2M 2 



AlTLO 



(AM 4 - uj 4 )L 2 



JU 2 



4M 2 ~ 



+ Vu7 2 



AM 2 ) 2 



(uj z + 2M l ) In 



AM 2 ) 

uj(uj 2 



3w 4 ) acosh 



+ (7M 4 + 2M 2 uj 2 



UJ 

2M 



AM 2 



M 3 



- (uj 2 + 6M 2 



+ 0(e,g 4 



(4.5) 
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where the function L 2 is defined as 



L 2 (x) = 4Li 2 (x) + 2Li 2 (-x) + [21n(l - x) + ln(l + x)]lnx . 



(4.6) 



The result in eq. (|4.5p agrees with the classic result from the literature pQ-[3]. The thermal 
correction, in turn, reads, 



8g 2 C A C F 
(Att) 3 uj 2 

9{uj)e(k 



AM 2 



U3 



2uj 



2u 2 k 2 < 



AM 2 



u{uj + 2k) 



+ (uj 2 + 2MVw(u + 2k)^uj(uj + 2k) - AM 2 
-2(w 4 - AM 4 + 2u)k{y? + 2M 2 ) + 2u?k 2 ^\ acosh 
,2 



uj(uj + 2fc) 
4M 2 



- 2M) f 



AM" 



2uj 



k 



2uj A k 



AM 2 



u)(u — 2k) 



+ (uj 2 + 2M 2 )y/u> (u - 2k)y/uj(u - 2k) - AM 2 



-2[uj a - AM 4 - 2ujk{oo 2 + 2M 2 ) + 2uj 2 k 2 ) acosh 
9[u - 2M) 



uj{uj — 2k) 
AM 2 



-2{uj 2 + 2M 2 ) ujyju 2 - AM 2 



+ .1. ( ^ 1 - I A / 1 + 2a.- 2 /.- 2 ) acosh ( 



0(e-? M ,g*) 



(4.7) 



where we have restricted to uj > (u < follows from antisymmetry, py(-uj) = — py(oj)). 
Eq. (|4.7p is our main result. 

A numerical evaluation of eq. (|4.7p , compared with the vacuum part in eq. (|4.5p , is shown in 
fig.[TJ We note that even though the thermal part is not exponentially suppressed for lo > 2M, 
it still only amounts to a small correction at phenomenologically interesting temperatures. 
On the other hand, the thermal part does possess the new qualitative feature that the result is 
non-zero below the threshold as well, where it is then the dominant effect; this can be traced 
back to reactions where a heavy quark and anti-quark annihilate into a gluon remaining inside 
the thermal medium, and a photon escaping from it. 

As an amusing remark, we note that while the next-to-leading order vacuum part is dis- 
continuous at the threshold, the next-to-leading order thermal part appears to be continuous. 
A similar pattern holds also for the scalar channel (fig. Ej): then the next-to-leading order 
vacuum part is continuous, while the next-to-leading order thermal part appears to have a 
continuous first derivative. These features are perhaps a manifestation of the fact that a non- 
zero temperature in general "smoothens" the spectral function; in a resummed framework, 
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Figure 1: The vacuum and thermal parts of the next-to-leading order correction in the vector channel, 
normalized by dividing with — ^uJ 1 g 2 C aG f ■ The vacuum part remains finite for uo — > oo (in units of 
the figure, its asymptotic value is 3/5127T 3 ), while the thermal part disappears fast for ui/M ^> 1. 

it may then not be surprising if any resonance peak of the vacuum result should disappear 
from the spectral function at high enough temperatures. 

To summarize, the characteristic feature of fig. [T] is a significant "threshold enhancement", 
due mostly to the vacuum part at T <C M. Within perturbation theory, this is to be 
interpreted as a first term of a series which, when summed to all orders, builds up possible 
quarkonium resonance peaks at uj < 2M. At the same time, the result of a resummed 
computation (to be discussed in more detail at the beginning of the next section) should 
extrapolate towards the perturbative one at some uj > 2M. 



5. Phenomenological implications 

We would now like to combine our result with that obtained within an NRQCD [312 EI] an d 
PNRQCD [32l [33] inspired resummed framework in ref. [§]. In order to do this, we need to 
pay attention to the correct normalization of the resummed result. In fact, the well-known 
(vacuum) normalization factor can be read off from eq. (14, 5ft : denoting 

Vuj 2 - 4M 2 

v = , (5.1) 
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M = 2 GeV, T = 400 MeV 



T 



T 



T 



resummed 
NLO QCD 
free theory 




0.0 - — 
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2.5 
co/M 
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Figure 2: A comparison of the near-threshold "resummed" result of ref. [S], matched to the "NLO 
QCD" expression of the present paper through an overall normalization factor, as discussed in the 
text. The difference of the NLO QCD result and the free theory result contains both the vacuum part 
and the thermal part; the magnitude of the latter is reflected by how much the curve deviates from 
zero below the threshold. 



the leading order vacuum expression can be expanded near the threshold as 



py(uS) 



lJ- 



6{uj - 2M) 



3C> 
8vr 



while the next-to-leading order vacuum result becomes 



Pv(uj) 



Ll)- 



8g 2 C A C F 6(LU - 2M) 



+ 0(v d 



Sv 



512tt 64vr 3 



+ 0(v 2 ) 



(5.2) 



(5.3) 



Since radiative corrections within a non-relativistic framework always contain a power of v, it 
is possible to account for the second term in eq. (|5.3p . equalling —g 2 CF/^ 2 times the leading 
term in eq. (15. 2j) . only by a multiplicative correction of the current |f| 



QCD 



g 2 c F 

2vr 2 



+ ... 



(5.4) 



In principle the coupling here should be evaluated at the scale ~ M [35], but in practice our 
resolution is low enough that we follow a simpler recipe (cf. next paragraph). In any case, 



3 The same relation is valid both for NRQCD and PNRQCD [34] 



11 



M = 2 GeV 



M = 4 GeV 



M = 6 GeV 



= 250 MeV 
= 300 MeV 
= 350 MeV 
= 400 MeV 
= 450 MeV 
= 500 MeV 



2.5 
(o/M 



= 250 MeV 
= 300 MeV 
= 350 MeV 
= 400 MeV 
= 450 MeV 
= 500McV - 



2.5 
(O/M 



2.2 
2.0 
1.8 
1.6 
1.4 
1.2 

• 1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



= 250 MeV 
= 300 MeV ~ 
= 350 MeV - 
= 400 MeV 
= 450 MeV ~ 
= 500 MeV - 



Figure 3: The phenomenologically assembled vector channel spectral function pv(w), in units of — oj 2 , 
for M = 2,4,6 GeV (from left to right). To the order considered, M is the heavy quark pole mass. 
Note that for better visibility, the axis ranges are different in the rightmost figure. 

the normalization factor is numerically significant, and its precise treatment plays a role; 
we actually do not impose it exactly, but rather search for a value minimizing the squared 
difference of the two results in the range (to — 2M)/M = 0.0 — 0.4, thereby also accounting 
for thermal corrections. This results in a normalization factor in the range 0.7 — 0.9, which 
indeed is the same ballpark as suggested by (the square of) eq. (|5.4| . given our choice of 
g 2 (cf. next paragraph). The "interpolated", or rather "assembled" result, is subsequently 

1 n l (assembled) / (QCD) (resummcd) \ A i r n .1 - i i i ■ l 

denned as p v = max(pj/ , Pv )■ An example for how the interpolation works 

in practice is shown in fig. [2J 

As far as the value of g 2 goes, no systematic choice is possible in the absence of NNLO com- 
putations at finite temperature. We follow here a purely phenomenological recipe, whereby 
g 2 is taken from another context where a sufficient level has been reached [36] , and take [37] 

^ 9^9.082^) ' *»iW,-3. (5.5) 

We also fix A^g ~ 300 MeV to be compatible with ref. [5]. It should be obvious that the 
subsequent results contain unknown uncertainties; still, the situation could in principle be 
systematically improved upon through higher order computations. 

The resulting full spectral function is shown in figs. [3l H] for various masses and tempera- 
tures, as a function of oj. The corresponding dilepton production rate from eq. (|2,2|) is shown 
in fig. Compared with the results in ref. [8] , the absolute magnitude of the rate has de- 
creased by about 10 - 30%, due to the inclusion of the normalization factor. We should again 
stress that particularly the charmonium case contains large uncertainties, and our results are 
to be trusted on the qualitative level only. 
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Figure 4: The phenomenologically assembled vector channel spectral function py(uj), in units of — oj 2 , 
for T = 250,350,450 MeV (from left to right). To the order considered, M is the heavy quark pole 
mass. Note that for better visibility, the axis ranges are different in the leftmost figure. 

6. Conclusions 

The purpose of this paper has been to compute the heavy quark contribution to the spectral 
function of the electromagnetic current at next-to-leading order in thermal QCD. The result 
consists of a well-known vacuum part, eq. (|4.5p . and a new thermal part, eq. (|4.7p . The 
thermal part is illustrated numerically in fig. [I] in comparison with the vacuum part. 

The thermal corrections in our result arise exclusively from the gluons with which the 
heavy quarks interact. Although these contributions are not exponentially suppressed, they 
turn out to be power-suppressed at large energies u ^> 2M: their general magnitude is 
0(g 2 T 2 ), and given that T < gM is the phenomenologically interesting temperature range 
(cf. eq. (|2.3p ). they can in principle be omitted in comparison with the next-to-leading order 
zero-temperature corrections, of 0(g 2 M 2 ). This also means that the asymptotic behaviour 
of the spectral function, needed as input for lattice studies, could (in the continuum limit) 
be extracted from the well-studied zero-temperature computations (see, e.g., ref [24]). At 
zero temperature the next-to-leading order correction could, perhaps, even be worked out at 
a finite lattice spacing. 

On the other hand, decreasing the energy towards the threshold, the thermal corrections 
become increasingly important. In fact, at next-to-leading order, the vacuum spectral func- 
tion vanishes at u) < 2M, while the thermal correction stays finite. The result emerges from 
phase space integrals associated with the energy constraint 5{u + k — E\ — E2), where uj is the 
photon energy; k is the gluon energy; and E\, Ei are the energies of a heavy quark and anti- 
quark. Graphically, the process corresponds to the annihilation of quarkonium into a gluon 
and a photon, the former of which remains within the thermal medium. Since large values of 
k are Boltzmann suppressed, the thermal corrections are substantial only for \co — 2M\ < T. 
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Figure 5: The physical dilcpton production rate, eq. (|2.2[) . from charmonium (left) and bottomonium 
(right), as a function of the energy, for various temperatures. The mass M corresponds to the pole 
mass, and is subject to uncertainties of several hundred MeV; we use the intervals 1.5. ..2.0 GeV and 
4. 5. ..5.0 GeV to illustrate the uncertainties. The low mass corresponds to the upper edge of each 
band. Compared with ref. [S], the main change is a 10 - 30% reduction of the overall magnitude. 

Combining our new results, valid far enough away from the threshold, with previously de- 
termined resummed expressions, valid close to the threshold, we have subsequently assembled 
phenomenological estimates for the spectral function in a macroscopic energy range (figs. [3l 
HJ. The corresponding dilepton production rate is shown in fig. Analogous results and 
plots for the spectral function in the scalar channel have been given in appendix C. The com- 
putations of the present paper play an important role in these plots particularly in that they 
fix the overall normalization of the assembled curves. We hope that these results can even- 
tually be incorporated in a simulation including an expanding and cooling thermal fireball, 
which would then allow for a direct comparison with the dilepton production rate measured 
in heavy ion collision experiments. 

We note, finally, that we have restricted to u > in this paper. There is a lot of interesting 
structure in the vector channel spectral function also around to ~ 0, related to the heavy quark 
diffusion coefficient. However, that structure is suppressed by exp(— (3M), and a non-trivial 
result also only arises at the order 0{a 2 s ) [38], so that our present computation at 0(a s ) 
cannot add anything to the known results |39j . 
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Appendix A. Intermediate steps for a master sum-integral 

We elaborate in this appendix on the steps outlined in sec. 13.31 The starting point is the 
expression in eq. f|3. 11 [) . 

A.l. Matsubara sums and the spectral function 

The first step is to carry out the Matsubara sums TJ2k n , ^Yl{ Pn }- The sum T^2 kn is 
complicated by the appearance of the functions II^, ITr in the gluon propagators. The reason 
for their introduction was that there could in principle be infrared divergences associated with 
the gluons; in the Euclidean formalism, these would come from small spatial momenta k for 
the Matsubara zero mode k n = 0, and could then be regulated by the fact that IIe(0, k) = 
> 0. Our strategy in the following will be to assume that there are no infrared divergences, 
whereby we can set TIt = He = 0; the absence of divergences will be verified a posteriori. 
Nevertheless, it has still been important to keep He ^ Ht in eq. (j3. IIP , because it could 
happen that the structure multiplied by \j{K 2 + 11^) — 1/(K 2 +TIe), which vanishes in the 
free limit, contains infrared sensitive parts which do not completely cancel against each other 
in the presence of He ^ Ht- 

In order to allow for an eventual introduction of He and 11^, we write for the moment the 
gluon propagators in the spectral representation, 



with k = |k|, the whole procedure is obviously just a rewriting of the decomposition 




(A.l) 



and carry out the Matsubara sum T ^2 kn with the kernel l/(k° — ik n ). In the free case, when 
the spectral function reads 



Piiee (k°, k) = — 5(k° - k) - 5(k° + k) 



(A.2) 



1 



1 



1 



1 



(A.3) 



2^2 j_ U2 



2/l k %k> 



+ 



k ~\~ ik' : 



'71 



■n 
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Note that the procedure is rather versatile and could also be interpreted as 



K 2 + U(K) 



U(K) 



K 2 + IL(K) 
dk° J > fre e(fc°,k) p(k°,-k) ( 
it \ k°-ik n 2k 



1 

K 2 



(A.4) 



1 



1 



+ 



1 



1 



1 



k - k° 



\ TXi \^ Lhjyi fh v t\j 72 rV li rn fi i\j brvy. 

where p(k°,k.) is the spectral function corresponding to — 11/ {K 2 + II), and we made use of 
the spectral function's antisymmetry in k° — > —k°. All the sums over k n are now with the 
same kernel as the one following from eq. (| A. 1 1) . The representation in eq. (|X4]) would be 
relevant for H^, in which case pe would have a pole at k° = k |25j . 

After this lengthy introduction, we are ready to carry out the sums. We describe the 
procedure in some detail for one of the master sum-integrals appearing in eq. f|3. 1 1 [) : for the 
others, the results are listed in appendix B. 
The case we choose to consider in detail is 

dk°^£ p(k°,k) 1 



S9(u}) = Disc 



Denoting 



vr tr K{p} k°-ik n A(P)A(P-Q)A(P-K) 



E p = VP 2 + M 2 , E p _ k = V(p-k) 2 + M 2 



(A.5) 
(A.6) 



we can rewrite the sums as 



{Pn} 



[*° - i*»K + E?H(p„ - u,„r + E|][(p„ - i„) 2 + f? J 



T 'EEEE 

k ri {p n } {r n } {Sn} 

P rP 



05,, 



r n —p 71 



n ,oPS. 



[fco - ik n ]\p 2 n + E 2 ][r 2 n + E 2 ][s 2 n + E 2 



dr/ dae 1 



o 



,Tr4 EEEE 



p— fcJ 



k° 



ik n p 2 +E 2 r 2 +E 2 s 2 +E 2 _ 



(A.7) 



kn { Pn } {r n } {s n } " ' -V ■ " ' -V ~n • ~p- k 

where in the last step we used a representation of the Kronecker delta-function, (38t n fi 
J^dTe itnT . The sums have factorized and can now be carried out: 

ra B (£; )e (<Tmod/3)fc , < amod/3 < (3 , (A. 



T Y1 k o 



±ir n r 

ri + E 2 



n F (E p 
2E„ 



= (/3-|rmod2/3|)£; _ \rmod2f3\E 



-(3 < t mod 2(3 < f3 , 



(A.9) 



where n F (u) = l/[exp(f3u>) + 1] and n B (u>) = l/[exp(/3u;) — 1]@ The subsequent integrals over 



t and a are elementary; we simply need to split f^da = f£ T da + fg_ T da, and note that in 



rP~r A 



rP 



4 The sum in eq. ()A.8|I is discontinuous at a = mod [3, and denning its value at the discontinuity re- 
quires care; although of no importance in the present context, we note that the expression with the correct 



antisymmetry in k° corresponds to the "average" , tib (k )(l + e< 3k )/2 = n B (fc°) 
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the latter range, |r + crmod2/3| = 2(3 — r — a. Setting e luJn @ = 1 after the integrations, the 
a> n -dependence of the result appears only in structures like l/(iu> n + ]TV Ei), and we can read 
off the discontinuity: 



Disc 



1 



(A.ll) 



= -TTd(u + J2 E i) ■ ( A ' 10 ) 

Implementing these steps in practice, and restricting the /c°-integral to positive values by 
making use of the antisymmetry of p(k°,h), we arrive at 

n Jk, P AE p E p _ k I 

8{u - 2E P ) - 8(u + 2E p ) (1 - 2n F1 ) x 

(A+\. + AI^)(1 + n B0 - n F2 ) - (All + A^)(n B0 + n F 
5(u - A ++ ) - 5(u + A ++ ) A++AI+ (1 + n B0 )(l - n F1 - n F2 ) + n F1 n F2 
5(u> — A ) — 5{u + A ) Al_A;7__ — n B0 (l — n m — n F2 ) + n F1 n F2 

ATiAli 



1 

2 El, 



5(uj - A + _) - <5(w + A + _) AliA^i n BO n F1 - (1 + n B0 )n F2 + n F1 n F 
5(oj - A_ + ) - <5(w + A_ + ) A++AI+ n B0 n F2 - (1 + n B0 )n F1 + n F1 n F 
To simplify the expression somewhat, we have introduced the shorthands 



A; + ovEL + 



p-k 



cr,r 



fir 



n F 



n F (£p) 



n F (E p . 



(A.12) 
(A.13) 



Note that the result in eq. (jA.lip is antisymmetric in u — > —uj, as must be the case. 

Inspecting eq. (jA.lip . we note the appearance of structures in the denominator, A^ etc, 

which look like they might vanish for some k, p. In fact, in one of the other master sum- 
integrals, even the structure l/{E p _ k — E p ) appears, which certainly vanishes, for 2pk = k 2 . 
It can be verified, however, that such poles always cancel between the various types of terms in 
the expression, and do not hinder the actual integration. (If integration variables are changed 
in a subset of the expression, p — > — p + k, to remove an apparent symmetry in E p <-> E p _/ t , 
then such terms do not in general cancel any more; nevertheless their contribution remains 
finite and correct if the poles are interpreted as principal values.) 



A. 2. Spatial momentum integrals 

The result so far, eq. (jA.lip . contains integrals with two types of delta-functions: ones with 
5(uj ± 2E p ), which we call "factorized" (fz) integrals, because the gluon momentum k does 
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not appear inside the 5-functions; and ones with more complicated 5-functions, which we call 
"phase space" (ps) integrals. In both cases, our strategy is to first carry out the integral 
over the quark momentum p = |p| and over the angle between p and k; the integral over 
the gluon momentum k = |k| is left for later (it is this integral which could potentially suffer 
from infrared divergences). 

We start by considering the phase space integrals, which are ultraviolet finite, so that we 
can set d = 3. In order to simplify the task, we ignore from now on terms suppressed by 
exp(— (3M) <C 1. This means that all appearances of n F {E p ) and n F (E p _ k ) can be omitted. 

Furthermore, restricting to uj > 0, we note that the delta-function 5(uj — A ) = 5(uj + 

E p + Ep_k — k°) can only be realized for k° > 2M, and will then lead to an exponentially 
small contribution due to the appearance of the Bose distribution n B (k°). The delta-function 
5(uj + A ++ ) = 5(uj + k° + E p + E p _ k ) does not get realized at all. Thereby only two of the 
eight delta- functions in eq. (jA.lip remain non-zero, and the integral simplifies to 

'dk° f d 3 k „Q f d 3 p 7T 



S 4 » 



(2vr)3 E p Ep_ k 

S(u -k°-E p - E p _ k )[l + n B (k°)]0(k°) + 6(u + k° - E p - E p _ k )n B (k°)0(-k )^ , (A.14) 
where 

' (A.15) 



4(A;0 + E p + E p _ k )(k0 -E p + E p _ k ) 



Fixing k and denoting z = —p-h/pk, so that E p _ k = y 7 p 2 + k 2 + 2pkz + M 2 , we can change 
integration variables from p, z to E p , E p _ k : 

d 3 p = 2vr / dpp 2 / dz = 2vr/ dE p / d£ p _ fc p p ~ k , (A.16) 
Jo J-i Jm J E~_ k * 



where E^_ k = \J p 2 ± 2pk + k 2 + M 2 . The hard task is to figure out when the (^-functions get 
realized. For 5 
which leads to 



p— k 

realized. For 5(u> — k° — E p — E p _ k ) this happens provided that E~ k < uj — k° — E p < E p _ k , 




k° < uj - 2M , k < vV - k ) 2 - 4M 2 , (A.17) 

4M 2 uj-k° k I AM 2 

- j- _ kQ)2 _ k2 <E p < + - y 1 - _ kQ)2 _ k2 . (A.18) 

In the case of the free gluon spectral function, with k° = k, these simplify to 

uj>2M , k < , (A.19) 



uj 




uj(uj — 2k) 



(A.20) 
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For 5{<jj + k° — E p — J5p_fc), we simply need to set k° — > —k° in eqs. (lATTll . (lA~l~8|) : putting 
subsequently k° = k, the explicit expressions read 

/ 4M 2 - u? \ 

uj>0, fc>max|0, - j, (A.21) 



oj + k k AM 2 „ oj + k k AM 2 fk , 

" ' - — ; — < E v < h-Wl : — - . (A. 22) 

u{u + 2k) p 2 2]j uj(uj + 2k) { ' 

Note also that the function <f> evaluates to — 1/[4w(uj — 2E P )] after integration over E p _k, for 
both delta functions in eq. (|A.14|) . 

Inserting the free gluon spectral function from eq. (|A.2p and using the simplified formulae 
from eqs. (jATT9]) - (|AT22|) . the integrals over E p _k and E p can now be carried out. For the 
thermal part, i.e. the one proportional to n B0 , this yields 



i^-4M 2 \ u(uj-2k) 



9(u - 2M)6[ k acosh 



2cu / V 4M 2 



+ 0(e" 



-/SAT*, 



The vacuum part, on the other hand, is given by the latter row of eq. (|A.23|) . but just without 
the function n B (k); then the final fe-integral is doable as well, and we end up with 



*4»ir = JJZTs^ -2M) 



ps 



(4vr 



(oj 2 -4M 2 )2 2M 2 -u 2 ( uj \ 

H Trt acosh — — 



4a; 2u 2 V 2M 7 



(A.24) 



Consider next the factorized integrals, i.e. the first term inside the curly brackets in 
eq. (|A.11|) . Again we start by integrating over p, z, and leave the integration over k for 
later. This time it is useful to view z as part of the k-integral, i.e. 

where d = 3— 2e. The factorized integrals are, in general, ultraviolet divergent, and necessitate 
keeping track of e ^ 0. As always, a helpful strategy is to add and subtract a simple infrared 
finite regulator, such as l/(fc 2 + M 2 ) a , where a is so chosen that the complicated expression 
becomes ultraviolet finite after the subtraction, and can be worked out at e = 0, while the 
ultraviolet divergent integral with the measure of eq. (|A.25P is taken over the simple regulator. 
In the complicated but ultraviolet finite integral, it is useful to change integration variables 
from z to Ep-k, using 



-i Ep_ k J E -_ h pk 



(A.26) 
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We should remark that in our particular example, S®, the trick of adding and subtracting 
a regulator is superfluous, given that the divergent integral can be directly identified as a 
known case, but in the general case we have found it to be very helpful. 

Now, because of the constraint 8(u — 2E p ) (for to > 0) in the factorized integrals, the 
integral over p can be carried out trivially. In fact, comparing eq. (jA.lip with (|B.3p . which 
defines a corresponding 1-loop integral (denoted by Si(u) and given explicitly in eq. (IB.4j) ). 
we arrive at 



'dA; 



7T 



2c 



d d k 



-p(fe°,k) 



(2TT) d 2E p _ k 

(Mi + A- 1 _)(n B0 + n f 



. (A.27) 



(A;^ + A_\){1 + n B0 -n F2 ) 

p=V^ 2 -4M 2 /2 

Let us first inspect the vacuum (T = 0) part hereof, i.e. the term without n B o or n F2 . Inserting 
the free gluon spectral function from eq. (lA.2p . the multiplier of Si(a>) becomes 



B = n 



2c 



d d k 



1 



1 



(2Tr) d 4kE p . k 
This can be compared with the integral 

d D K 



+ 



k + E p + Ep_k k — E p + Ep_k 



p=Vw 2 -4M 2 /2 



B (P z ;0,M z ) 



2c 



(1 



(2tt) d K 2 [{P- K) 2 + M 2 } 



:'2c 



d d k 



1 



ipo + k + E. 



+ 



1 



p—k 



-ipo + k + E, 



p-k 



(A.28) 

(A.29) 
, (A.30) 



(2vr) d AkE p _ k 

where we denoted K = (&o,k) and carried out the integral over ko. In other words, Bo 
B (P 2 ;0,M 2 ), where 



Pe 



-iE p ,pe) 



p=VaJ 2 -4M 2 /2 ' 



1 A 



-M 



and e is a unit vector; the value of this standard vacuum integral reads 



B (-M 2 ;0,M 2 ) 



1 



(4vr) 



(A.31) 



(A.32) 



Combining this with eq. (|B.4|) . the factorized vacuum part becomes 
S»|™ c = S 1 (u;)Bo(-M 2 ;0,M 2 ) 



0(w - 2M)- — - S J tanh' 



4u;(47r) 



1+1 » 

e n M 2 (w 2 -4M 2 ) 



+ 4 + 0(e) 



(A.33) 



For completeness, we have even kept exponentially small thermal terms in the coefficient of 
1/e, given that it is useful to crosscheck the exact cancellation of ultraviolet poles; after this 
check, we set tanh(/3w/4) = 1 + 0(exp(-/3M)), given that u > 2M. 
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Consider then the thermal part of eq. (|A.27|) . Again, we omit exponentially small terms 
~ exp(— (3M), and use the free gluon spectral function. Because of the remaining factor n B0 , 
the &:-integral is exponentially convergent, and we can set e = 0. Employing eq. (|A.26P the 
thermal part becomes 



s 4 V) 



{At:) 2 
1 



dkkn B (k) 
1 



dz 



E, 



p-k 



+ 



1 



k + E n + E p _ k k-E ri + E, 



p—k 



k-E n - E, 



p-k 



k + E n -E. 



p-k 



p=Vu 2 -4M 2 /2 



(4vr) 2 p J 



dkn B (k) x 



x In 



(k + E p + E+_ h )(k -E p + E;_ k ){k -E p - E+_ k ){k + E p - E+_ k ) 



{k + E p + E;__ k )(k -E p + E p _ k )(k -E p - E;_ k )(k + E P - E^ 



P=Vuj 2 -4:M 2 /2 

(A.34) 



Making use of 



(k + aE p ) 2 - {E T p _ k f = 2k[oE p - rp] , a,T = ± 



it can be seen that the argument of the logarithm evaluates to unity. Hence, S®(uj) 
To summarize, combining eqs. (|A.23|h (|A.24p . (|A.33p . we get 

s°M = s2(u)\ T ps + s2(u)\ v ;; + 5 4 V 



vac 
fz 



(A.35) 



0. 



(A.36) 



The other master sum-integrals can be worked out in the same way, and the final results are 
listed in appendix B. 
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Appendix B. General results for all master sum-integrals 



We collect in this appendix the results for all the master sum-integrals entering the computa- 
tion, obtained with the methods explained in appendix A. In each case, we list the definition 
of the sum-integral; an intermediate result obtained after carrying our the Matsubara sums 
and taking the discontinuity; and the final result, obtained after restricting to the free gluon 
spectral function, omitting terms suppressed by exp(— (5M) (except from the ultraviolet di- 
vergences), and carrying out the final spatial integrations. As before, the integration measure 
for the spatial integrations is defined as 



= fj, 



2c 



d 3 " 2e p 



(27T) 



3-2e 



(B.l) 



and p, 2 = 4nfi 2 e~^ E denotes the MS scale parameter. To simplify the expressions somewhat, 
we also make use of the shorthands listed in eqs. (|A.12p . (|A.13|) . The subscripts "ps" and 
"fz" denote "phase space" and "factorized" integrations, respectively, in the sense of sec. IA.2I 



B.l. Si 

The sum-integral Si is defined as 

Si(u) = Disc 



{P} 



A(P)A(P-Q) 



Q={-iui,0) 



Carrying out the Matsubara sum and taking the discontinuity leads to 



7T 

AE 2 
p v 



l-2n F (£p) S(uj - 2E P ) - 8{uj + 2E p 



(B.2) 



(B.3) 



The remaining integral is trivial due to the 5-function and, restricting to u > 0, we arrive at 



5i(w) =9(u-2M) 



(uj 2 - 4M 



tanhf^^ 



1 + e In 



2 -AM 2 



+ 2 + 0(e 2 



(B.4) 



B.2. S 2 

The sum-integral S2 is defined as 
52 (to) = Disc 



{P} 



A 2 (P)A(P-Q) 



=(—iu,0) 



It is easy to see that ^2 = -\&Si/dM 2 . Therefore, from eq. (lR4l) . we obtain 



S 2 (cu) = 6{uj- 2M) 



(uj 2 - 4M 2 
16ttuj 



„[* 



1 + e In ■ 



uJ 



2 -AM 2 



+ 0(e 2 ) 



(B.5) 



(B.6) 
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B.3. S 3 

The sum-integral S3 is denned as 



Sx(uj) = Disc 



dk°j: p{k°,k) 1 

~T K{p} k° - ik n A(P)A(P — Q - K) 



Q=(-iu,Q) 



(B.7) 



Performing the Matsubara sums, taking the discontinuity, and making use of the antisym- 
metry of p(k , k) yields 



Jo 71 Jk, P Ah p b p - k y 



p -*j-lpj-lp—k 

8{uj - A ++ ) - 5(ijj + A ++ ) (1 + n B0 )(l - n F1 - n F2 ) + n F1 n F2 

5{uj — A ) — 5(lo + A ) — n B0 (l — n F1 — n F2 ) + n F1 n 

6(u) - A + _) - 5(uj + A + _) n BO n F1 - (1 + n B o)ra F 2 + n F1 n F2 

£(cj - A_ + ) - 5(u + A_ + ) n BO n F2 - (1 + n B0 )n F1 + n F1 n F2 



(Bi 



+ 



Inserting the free gluon spectral function, and omitting exponentially small terms, yields 



(47T) 



0(w - 2M) 



(aj 2 -4M 2 )5 ( W 2 + 2M 2 ) M 2 (M 2 -^ 2 ) 



8lj 



+ 



+ / dkkn B (k) 
Jo 



e { .)0(k- iM2 -^), 1- / M2 M 

1 2uj JV uj(uj + 2k) 



roj 2 - 4M 2 
+0(a,-2M)fl( - 2 J -A-),/ i 



4M 2 



— 2fe) 



(B.9) 



B.4. 5° 



The sum- integral is defined in eq. (|A.5P ; its value after the Matsubara sums is given in 
eq. (|A.11[) : the result after the phase space integrals is the sum of eqs. (|A.23p . (|A.24|) . (|A.33|) . 



B.5. Si 



The sum-integral S4 is defined as 



SI(uj) = Disc 



7T f K{p} k° - ik n A(P) A(P - Q) A(P - tf) 



Q=(-icj,0) 



(B.10) 
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Performing the Matsubara sums, taking the discontinuity, and making use of the antisym- 
metry of p(k°,k) yields 



^ r dk ° f (ifi ^ wk0uj / 

JO * Jk,p AE p E p - k [ 



(B.ll) 



2 E n 



6{u - 2E P ) + 8{u + 2E p ) (1 - 2n F1 ) x 



x 
+ 
+ 
+ 

+ 



- Ali)(l + n B0 - n F2 ) + (All - A7l)(n B0 + n F2 ) 



- A ++ ) + S(u + A ++ ) 

5{uj - A ) +5{uj + A__) 

S(oj- A+_) + A+_) 

<y(w - A_+) + <5(u; + A_+) 



ATiAli 



AllA;l 
AllA+l 



A'lAzl 



(1 + n B0 )(l - n F1 - n F2 ) + n F1 n F2 

— ^Bo(l ^Fl ~i~ ^P1 ? T'F2 

(1 + n B0 )n F2 + 7T- F l/Z F2 

^"B0^F2 (1 ~\~ ^Bo)^Fl ~\~ npiTlp2 



Inserting the free gluon spectral function, the ultraviolet divergent factorized vacuum part 
reads 



vac 
fz 



-6{uj - 2M) 



UJ{UJ 



2 -4M 2 )! j(3uj- 



16(4tt) 2 



■ tanh 



(?) 



1 , A 4 

e n M 2 H -4M 2 ) 



+ 3 + 0(e) , 
(B.12) 



where in the coefficient of the divergence we have accounted even for exponentially small 
terms. The vacuum part from the phase space integrals reads 



9{u> - 2M) 
(4vr) 3 



3{u 2 -AM 2 ) 2 (2M 2 -J 1 
32^ 



+ 



4 - 4a; 2 M 2 + 6M 4 



8lo 2 



acosh ( ) 



while the thermal parts amount to 
1 



si hi: 



( 4?r )3 



d/c/cn B (£;)6>(cj - 2M) 



-2 acosh 



2M 



+ 0(e 



■*( W -2M)*( ^ 



*H^- 4MW u OS h^ +2fc) 



2uj 

kj acosh 



4M 2 



\2mj\ ' 

(B.13) 



(B.14) 
(B.15) 



uj 2 - AM 2 ,\ , luj{uj-2k) 

4M 2 



+ C»(e-^ M ) 
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B.6. S% 



The sum-integral «Sj is denned as 



5*4(0;) = Disc 



K{P} 



K 2 A(P) A(P - Q)A(P - if) 



=(— iu,0) 



(B.16) 



Because of the ultraviolet divergent factor in the numerator, the use of the spectral repre- 
sentation requires care in this case, and we rather proceed directly with the sum, having 
gone over into free gluons to start with. Carrying out the shift K — > P — K , the summation 
f actor izes, 



Disc 



1 



{P} 



A(P)A(P-Q) 



=(-tw,0) 



1 



{K} 



A(K) 



Si(u;)/o(M 2 ) , (B.17) 



where Si(u) is given in eq. (|B.4|) . while Iq, is a basic tadpole integral generalized to finite 
temperature. In fact, the finite temperature effects in Jo are exponentially small and can be 
omitted: 



I (M 2 







M 2 


L 2fc 


1 - 2n F (E k ) 


(4vr) 2 



7^ + ' 



+ <D(e,e~P M ). (B.18) 



Keeping exponentially small terms in the coefficient of the divergence, though, we arrive at 



S 2 (u)) = -0(u-2M) 



(u 2 -4M 2 )2M 2 



4o;(47r) 



tanh( ^— 
4 



1 , A* 4 

e U M 2 (co 2 -4M 2 ) 



+ 3 + 0{e,e- l3M ) 

(B.19) 



B.7. S% 



The sum-integral S9 is defined as 



5g(w) = Disc 



dk°j: p(k°,k) 



(B.20) 



=(-iw,0) 
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Performing the Matsubara sums, taking the discontinuity, and making use of the antisym- 
metry of p(k ,k) yields 

coo jlO 



2E. n 



5{oj - 2E P ) - 5{uj + 2E P ) 



:i-2n P1 ) x 



(a;^ + a:^)(a;^ - a:^ + e^){i + n B0 - n F2 ) 
+(a:^ + a;1)(a:1 - a;^ - E^)(n w + n F2 ) 



2E. n 



5(uj — 2E P ) — 5(uj + 2E P ) (1 - n F1 )n P1 x 



(A + ^ + A_*_)(l + n B0 - n P2 ) - (A_i + A + i)(n B0 + n P2 ) 



+ 



+ 



2 E n 



5' {uj - 2E P ) + 5\uj + 2E p ) 



[I - 2n F1 ) x 



+ 



+ 



+ 



(A + ^ + AI^)(1 + n B0 - n P2 ) - (All. + A;i)(n B0 + n P2 ) 

s(u - a ++ ) - s(u + a ++ )1 a^a:^(a:^ - a;^) x 



(1 + n B0 )(l - n F1 - n F2 ) + n F1 n F2 



5(lo - A__) - 5{uj + A__) AliAl* (Al* - All) x 



-n B0 (l - n F1 - n P2 ) + n F1 n F 



<5(w - A+_) - <5(kj + A+_) Al^AT^A^ - All) 



^BO^-Fl 

(1 + n B0 )n P2 + 



5(cc - A_ + ) - 5(u + A_+)j a^a:^(a:^ - a-\j 



(B.21) 



^B0^-F2 (1 "I - ^Bo)^Fl ~i~ W F i7ip 



In the factorized part, it is useful to carry out a partial integration in order to remove the 
structure 5'(u - 2E p ) + 5'(u + 2E P ): 



I 



(B.22) 



-/ 



d d p 

J2^y 



5{uj - 2E p 



(d-2)E p g | Epdg | 1 d(E p g) | (p + kz)E p 8g 



2p 2 



2p dp 2E P dEp 2pEp_ k dE p _ k J 

The subsequent steps proceed as described in appendix A. 

In contrast to S^, S4, however, it is not possible to give separate closed expressions for 
the factorized and phase space vacuum parts of S®, because the integrals are logarithmi- 
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cally divergent at the lower limit of the /c-integration (in the thermal case, they are linearly 
divergent). Yet the sum is finite, and inserting the free gluon spectral function, we get 



+ 



4a;(4vr) 3 (u; 2 -4M 2 )2 I V 4 
w 2 -4Af 2 1 w(w 2 -4M 2 ) (w 2 -4M 2 )5(4M 2 -3w 



e ' '" .M-l - 1.1 f-M 



+ 2 



(B.23) 



M 2 



In 



M 3 



+ 



wM 2 



acosh( — )|>+0(6), 



where in the coefficient of the divergence we have accounted even for exponentially small 
thermal corrections. For the thermal part proper we obtain 



Sg(")f 



1 



4w 2 M 2 (4tt)- j 



(B.24) 



0(u;)0(A;- — — — ) + 2k ) V^iu + 2fc) - 4M 2 



2uj 

+6{u- 2M)e( - ~ 4M - k) V^(w - 2k)y/uj(u - 2k) - 4M 2 



2cj 



-6{u> - 2M) x 2w\/w 2 - 4M 2 



The last line, which originates from the factorized integrals, subtracts the values of the first 
two lines at k = (for u > 2M), rendering the integral infrared finite. 



B.8. S% 



The sum-integral S® is defined as 



S£(u) = Disc 



dk°j: p(k°,k) 



(P • k) 5 



n T K{p} k°-ik n A 2 (P)A(P-Q)A(P-K) 



=(-iw,0) 



(B.25) 



Carrying out the Matsubara sums proceeds precisely like for S*P, and leads to an expression 



like eq. (|B.21h ; it is also again useful to carry out the partial integration in eq. ()B.22j) . The 
subsequent steps lead to the vacuum part 



0(w - 2M){u? -4M 2 ): 
8w(4vr) 3 



tanh 



/3u 



1 , A 4 

e n M 2 (w 2 - 4M 2 ) 



+ 1 



„ , cu(u; 2 - 4M 2 ) 2(7u 2 - 8M 2 ) , / u 
41n „„ + -i > acoshf — 



M 3 



WW 



4M 2 ); 



+ 



2w 



(w 2 -4M 2 )2 



(B.26) 
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Here the function 

a(v) = 



00 dx 
x 



9(v 2 -x)(l - x)ln 



\/l — x + Vv 2 — X 



VI — x — \f' 



— X 



(1 + x + y/1 + 2vx + x 2 )(-l + x + Vl - 2ra + x 2 )l 

t In - , (rJ.zY) 

(1 + x + VI - 2vx + x 2 )(-l + x + VI + 2t>x + x 2 ). 

where the integration variable x is related to k through k = xuj/2, is finite, but we have not 
bothered to work out its analytic expression, given that it does not appear in our final result. 
The thermal part reads 



1 



2u; 2 (4vr)~ Jo 



k 



(B.28) 



-^u{uj + 2k)^u{u + 2k) - AM 2 + u(u + 2fc) acosh 



+ 0(w-2M)0 



fc x 



Vw(w - 2k)y / uj(uj - 2k) - AM 2 + uj{uj - 2k) acosh 



— 2k) 
AM 2 



+ 6{u-2M)x 2uj^/uj 2 - AM 2 - 2u? acosh (j^j } + 0[e~ pM ) . 

The last line, which originates from the factorized integrals, subtracts the values of the first 
two lines at k = (for u> > 2M), rendering the integral infrared finite. 

B.9. Si 



The sum-integral S 2 is defined as 

S 2 (u) ee Disc ^ 1 



K 2 



K{P} 



K 2 A 2 (P)A(P - Q)A(P - K) 



Q=(-ia;,0) 



(B.29) 



Because of the ultraviolet divergent factor in the numerator, the use of the spectral represen- 
tation requires care in this case, and we rather proceed directly with the sum, as in the case 
of S 2 . Carrying out the shift K —> P — K ', the summation factorizes, 



S 2 (u>) = Disc 



1 



{P} 



A 2 (P)A(P-Q) 



Q=(-iu;,0) 



f, ■ X ^- ) =S 2 (u J )Io(M 2 ), (B.30) 
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where S^w) is given in eq. (|B.6p . while Iq is given in eq. (|B.18p , Keeping exponentially small 
terms in the coefficient of the divergence, we arrive at 



4o;(c; 2 -4M 2 )2(47r) 3 



1 



IT 



e +ln M 2 (w 2 -4M 2 ) 



+ l + 0(e,e- f3M ) 



(B.31) 



B.10. 5{? 

The sum-integral <Sg is defined as 



Sg(cj) = Disc 



dfc° 



P(fc°,k) 



=(— iw,0) 
(B.32) 

Performing the Matsubara sums, taking the discontinuity, and making use of the antisym- 
metry of p(k , k) yields 



(B.33) 



8£ P 2 



- 2£ p ) - 6(u + 2£ p ) (1 - 2ra P1 ) x 



ATi 



+ 



A 



-i 



E n — E„ 



(n B0 + n E 



p J-'p—k 

+ — ~ )(1+ n B0 — n F2 ) 



+ 



E p + E„_k E„ — Ep_k 



5(uj - 2£p_ fe ) - 8(u + 2E p _ k ) (1 - 2n F2 ) x 



A 



-i 



+ 



A 



-i 



Ep-k + E p E p _ k — E, 
All 



a: 



+ 



E r ,-k + En E, 



E r 



(n B0 + n E 



(1 + n B0 - n F1 ) 



+ 



5{uj - A ++ ) - 8{u + A++) 
6(u - A__) + A__) 

5(w - A + _) - 5{lo + A + _) 

- A_ + ) - 5{uj + A_ + ) 



ATiAli Ali 



a: 2 a;^a:^ 

AT 2 At! Ali 



AI 2 ATIAIT 



(1 + n B0 )(l - n F1 - n F2 ) + n F1 n F2 

^Bo(l ^Fl ^F2j 4" ^F1^F2 

n B0 "-Fi - (1 + n B0 )n F2 + n F1 n, 
« B0 n F2 - (1 + n B0 )n F1 + n F1 n, 



In the factorized part, a change of integration variables p — > k — p allows trivially to change 
the structure with <J(cj — 2E p _j t ) — 5(uj + 2E p _jS) into the familiar one with 6(u — 2E P ) — 
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5(lo + 2E P ). (The only complication is that then the difference 1/(E P — E p _k) needs to be 
interpreted as a principal value.) The subsequent steps proceed as described in appendix A. 

Like with S®, it is again not possible to give separate closed expressions for the factorized 
and phase space vacuum parts of S®, because the integrals are logarithmically divergent at 
the lower limit of the ^-integration (in the thermal case, they are linearly divergent). The 
sum is infrared finite, however, and inserting the free gluon spectral function, yields [3] 



0{uj - 2M) 
w 2 (4vr) 3 2 



AM 2 \ 



to 



AM 2 



(B.34) 



where the function L 2 is defined in eq. (14, 6p . For the thermal parts we obtain, omitting 
exponentially small terms, 

n B (k) 



w 2 (4tt) 



dfc- 



A: 



6{uj)6[k 



4M 2 



to 



2to 



acosh 



u(ui + 2k) 
AM 2 



-0(u - 2M) 9 



CO 



2 -AM 2 



2uj 



k ) acosh 



oj(u> — 2k) 



AM 2 



-0(u - 2M) x 2 acosh 



uO 



+ 0(e 



(B.35) 



2M J 

The last line, which originates from the factorized integrals, subtracts the values of the first 
two lines at k = (for 00 > 2M), rendering the integral infrared finite. 



B.ll. S% 

The sum-integral S* is defined as 



S%(u) = Disc 



dA; ^ p{k°,k) 



(p-k) 



vr J K{p} k ° ~ ik n A(P)A(P - Q)A(P - K)A(P - Q - K) 



Q=(-iuj,0) 

(B.36) 

Carrying out the Matsubara sums proceeds precisely like for S$, and leads to an expression 
like eq. (|B.33p . The subsequent steps lead to the vacuum part 

s 6 V) 



(lo - 2M) 



to 2 (Air) 



-M 2 L 2 



10 — \fbJ 
LO + V 00 



AM 2 \ ... , 
H 8 

1-AM 2 ' 2 H 



Vuj 2 - AM 2 



00 



+ 0(6). 

(B.37) 



Here the function 
r°°dx 
Jo x 



(5{v) 



+ 



6{v 2 - 



x)y/l — x\f v 2 — x — v + - + %vx + x 2 — yl — 2vx + x 2 



In 



(1 + \/l + 2vx + x 2 )(-l + s/l - 2vx + x 2 ) 



(1 + \/l-2vx + x 2 )(-l + y/1 + 2vx + x 2 ) 



(B.38) 
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where the integration variable x is related to k through k = xu/2, is finite, but we have not 
bothered to work out its analytic expression, given that it does not appear in our final result. 
The thermal part reads 



i 



2o; 2 (4tt) 



"«-<*> 



k 



(B.39) 



^oj(uj + 2k)y/u(u + 2k) - AM 2 - 4M 2 acoshW 

/w 2 -4M 2 

+ 0(u;-2M)0 fe | x 

\ 2uj 



y/uj{vo - 2k)y / oj(uj - 2k) - AM 2 - AM 2 acosh 
+ 0(w - 2M) x 



4M 2 



2w\/^ 2 -4M 2 + 8Af 2 acosh (^^ 

The last line, which originates from the factorized integrals, subtracts the values of the first 
two lines at k = (for u > 2M), rendering the integral infrared finite. 



B.12. S 2 

The sum-integral S 2 is defined as 

1 K 2 



Sq(uj) = Disc 



K{P} 



K 2 A(P)A(P - Q)A(P - K)A(P - Q - K) 



(B.40) 



Like with 5| and S 2 , we proceed directly with free gluons rather than using the spectral 
representation. Carrying out the shift K — > P — K , the summation factorizes, 



Sg(u) = Disc 



2 Re 



1 



{P} 



A(P)A(P-Q) 



=(-iw,0) 



1 



A(K)A(K - Q) 



=(-iw,0) 



{^} 



1 

A(P)A(P-Q) 



x Disc 



=(— iw,o) 



A(K)A(# - Q) 



=(— iw,0) 
(B.41) 



where Re[...] denotes the regular (non-discontinuous) part, while the discontinuous part can be 
identified with the function S\{u), given in eq. (|B.4|) . The Matsubara sum in the regular part 
can be carried out as before; the only difference with respect to the procedure in appendix A 
is that taking the regular part after the substitution uj n — ► — iuj yields a principle value rather 



31 



than a delta-function: 
Re 



{P} 



A(P)A(P-Q) 



1 



4P| 



P 



1 



u + 2E r 



- P 



1 



u - 2E r 



1 - 2n F (E p 



(B.42) 



It is seen that the finite-temperature effects continue to be exponentially suppressed. The 
zero-temperature part, on the other hand, equals the real part of the general function Bq, 
another special case of which was met in eq. (|A.29|) : 



Re 



"J 



d D P 



(2tt) d (P 2 + M 2 )[(P - Q) 2 + M 2 } 



Re 



=(— iu;,0) 



B (V;M 2 ,M 2 ) 



(Air) 



1 fl 2 2(w 2 -4M 2 )2 

— h In — — r + 2 acosh 

e M z uj 



2M 



+ 0(e) 



, w>2M. (B.43) 



Combining this with Si(u>), and keeping the exponentially small terms in the coefficient of 
the divergence, we arrive at 



(B.44) 



1 + 1 ^ 

e n M 2 (w 2 -4M 2 ) 



+ 4 



2(w 2 -4Af 2 ); 



arosh|^)+0( e ,e-^^ 



B.13. 52 



The sum-integral S 2 is defined as 



Sl{uj) 



Disc 



K 2 [p 2 - (p • k) 2 ] 



K{P} 



K 2 A(P)A(P - Q)A(P - if)A(P -Q-K) 



(B.45) 



The summation factorizes into two independent parts like for S 2 ; however, the spatial in- 
tegrations do not factorize due to the additional structure in the numerator. Therefore the 
evaluation is somewhat more involved, yet the general techniques introduced in appendix A 
yield a solution: 



1 , a 4 

e n M 2 (uj 2 - AM 2 ) 



+ 4 + 



x 

J 2 + 2M 2 
3(lo 2 - AM 2 ) 



(B.46) 



2(w 4 - 6w 2 M 2 + 12M 4 ) 
uj{uj 2 -AM 2 )1 



UJ 



M ^h[^j+0(e,e^ M ) 
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B.14. Si 



The sum-integral Sq is denned as 



SZ(u) = Disc 



tJ 



(K 



2\2 



K{P} 



K 2 A(P)A(P - Q)A(P - K)A(P - Q - K) 



Q=(-iu;,0) 



(B.47) 



Like with S 2 , <S 2 and Sg, we proceed directly with free gluons rather than using the spectral 
representation. Cancelling one K 2 and carrying out the shift K — > P — K, we get 



S 4 (w) = Disc 



A(P) + A(iT) - 2(M 2 + P ■ K) ' 
{Kp} A(P)A(P-Q)A(K)A(K-Q)\ Q=( _ (J , 



Another change of integration variables, P — > Q — P, shows that 



fp} A(P) A(P - Q) 2 ^ {p} A(P) A(P - Q) 
and similarly for the term with Thereby we arrive at 



(B.48) 



(B.49) 



5 6 4 H = 2S 4 2 H + -(^ 2 -4M 2 )S 6 2 H 



0(w - 2M) 



(to 2 - 4M 2 



4w(47r) 



■tanh^^^ x 



(w 2 - 6M 2 ) 
+ (kj 2 - 4M 2 ) 



I jU 4 

e + n M 2 (w 2 -4M 2 ) 



+ 3 



1 - 



2(lo 2 - 4M 2 



acosh 



a; 



— 

2M 



+ 0(e,e 



(B.50) 



(B.51) 
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Appendix C. Spectral function in the scalar channel 



For completeness, we have worked out the spectral function corresponding to the scalar 
channel with the same methods as described above for the vector channel. It seems, though, 
that the physical significance is not clear in the scalar case: no direct relation to an observable, 
in the sense of eq. (j2.2|) . has been worked out, as far as we know, and the computation as 
such appears to possess a number of ambiguities. In particular, the scalar density requires 
renormalization, and the renormalization factor cannot be uniquely specified; moreover the 
resummation of the spectral function within a potential model near the threshold appears to 
lead to ambiguities [8]. Nevertheless, on the lattice the correlator of (bare) scalar densities 
can be treated on the same footing as that of the vector currents [211 122] . 

Concerning the first of the issues, namely renormalization, the method we choose is to 
consider the object 

S = M<P$j>, (C.l) 

where is essentially the bare quark mass defined in eq. (|3.7p . only with a possible 

additional constant as a "probe", 

r(*A 2 _ ,,2 ^C F M^(l . ... P 2 . 4 



We then define 

ps{u) = f°° dteiUlt J d3 " 2ex (I [5(t,x),5(0,0)]\ , (C.3) 

which turns out to be finite. Starting again at 1-loop level, and omitting Q-independent 
terms which are killed by the discontinuity in eq. (|2.5p . we get 

<Z> - » " ""•*■] " ' < C ' 4 > 

The counterterm graph yields 

r ^ . , . 12g 2 C A C F M 2 f (\ R 2 4\ 

= [Q - mdep - ]+ { A,) 2 {U +ln ^ + 3j (C - 5) 
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Finally, the genuine 2-loop graphs add up to 
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(C.6) 



Again any dependence on the gauge parameter £ has disappeared, and P^ v of eq. (|3.2|) could 
have been replaced with 5^ u — P^ u . 

We note that the master sum-integrals appearing in eq. (|C.6|) are a subset of those in 
eq. (|3.1ip . Therefore the discussion in sec. 13.41 continues to hold, and there are no infrared 
divergences in the result, so that we can set TIt = n# = in eq. (|C.6p . The full result can 
now be written as 
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(C.7) 



We have set here e — > whenever the master sum-integral that it multiplies is finite. 

Now, the structure of eq. (|C.7p reveals an ambiguity with regard to the treatment of the 
"resummation" of thermal mass corrections, which in the vector case lead to eq. (|4.3p . In the 
vector case, the need to resum is unambiguous, because anything else than eq. (|4.3p would 
lead to a thermal correction diverging at the threshold. In the scalar case, we do not have 
this guidance: terms multiplied by T 2 vanish as 6{uj — 2M)(w — 2M)s at the threshold, 
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being thus subdominant with respect to the leading thermal corrections which remain non- 
zero. Nevertheless, we would like to apply a "universal" thermal resummation, i.e. precisely 
eq. (|4.3p ; however, it may be questioned whether it is valid to do this also in the term M 2 , 
coming from the ("ultraviolet-completed") definition of the scalar current, or only in more 
infrared sensitive parts. It seems to us that this question can be fully settled only through a 
next-to-next-to-leading order computation; in the following, we assume that the resummation 
of eq. (14.31) is only carried out in the Lagrangian, not in the definition of the scalar density. 
If so, a redefinition of the mass according to eq. (14.3P leads to the modified result 



p s (u) = 2C a M 2 (uj 2 - 4M 2 )Si(w) + Ag 2 C A C F M 2 
T 2 , x 3<5 
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+ A(uj 2 - AM 2 )Sl(uj) + 4(1 - e)S\(uj) - 2(uj 2 - AM 2 ) \2M 2 S%(uj) - (1 - e)Sf(w) 
+ (uj 2 - 2M 2 ){u? - 4M 2 )S 6 °(w) - \euj 2 + 4(1 - e)M 2 } S%(u)\ + 0(e) . (C. 



Unfortunately, the issue of what is resummed is not insignificant in the sense that the differ- 
ence between eqs. (|C.7P and (|C8[) is numerically of 0(1) for uj > 2M. 

Inserting the explicit expressions for the functions Sf (uj) from appendix B into eq. (|C.8p . 
the final result for the vacuum part reads 
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where the function L 2 is defined in eq. (|4.6|) . Let us note that although similar to the 
vector channel spectral function in eq. (|4.5|) at first sight, eq. (|C.9p has also some significant 
differences; in particular, logarithms of ui/M do not cancel at u 3> M any more, but the 
asymptotic behaviour becomes 



PS Hp 



u,»M 3g 2 C A C F uj 2 M 2 
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(C.IO) 



As the dependence on 5 and <5's definition through eq. ()C.2p show, the logarithm is in some 
sense related to the need to renormalize the scalar density and its correlators. 
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Figure 6: The vacuum and thermal parts of the next-to-leading order correction in the scalar channel, 
normalized by dividing with Alo 2 M 2 g 2 C aC f ■ The vacuum part can become negative because the 
bare scalar correlator is multiplied with a bare mass parameter, cf. eqs. (|C.lj) . (IC.2|) ; the constant 5 
illustrates how strong the dependence on the renormalization convention is. 



The thermal correction, in turn, reads 
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where we represented T 2 as 7r 2 T 2 = 6 / °° dkkn B (k). 
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M = 2 GeV 



M = 4 GeV 



M = 6 GeV 




Figure 7: The phenomenologically assembled scalar channel spectral function ps(lo) : in units of 
w 2 M 2 , for M = 2,4,6 GeV (from left to right). To the order considered, M is the heavy quark pole 
mass. Note that for better visibility, the axis ranges are different in the rightmost figure. As discussed 
after eq. (|C.14|) . we are not confident that these plots have a definite physical significance; the figures 
are meant for illustration only. 

A numerical evaluation of this result, compared with the vacuum part of eq. ([Cl9]h is 
shown in fig. [6J For small oj the thermal part appears to be somewhat more significant 
than in the case of the vector channel; this is because there is a cancellation of positive and 
negative contributions in the vacuum part, before the negative terms take over at large u (cf. 
eq. (jC.lOP ). The thermal part, in contrast, stays positive and vanishes rapidly at large uj. 

We wish to draw attention to the amusing feature, already mentioned at the end of sec. [U 
that while the next-to- leading order vacuum part is continuous, the next-to-leading order 
thermal part appears even to have a continuous first derivative. In the vector channel, in 
contrast, the next-to-leading order vacuum part is discontinuous at the threshold, while the 
next-to-leading order thermal part appears to be continuous (cf. fig. [T|). In other words, the 
thermal part seems always to be one degree smoother than the vacuum part. 

In order to now combine our result with that obtained within a resummed framework in 
ref. [8], we need to match the normalizations, in analogy with eq. (|5.4p . Indeed, employing 
the notation of eq. (|5.ip . the leading order vacuum result in eq. (|C,9|) becomes 



3 



oj 2 M 2 

while the next-to-leading order result can be expanded as 
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(C.13) 



Since radiative corrections within a non-relativistic potential model always involve a power of 
v, it is possible to account for the second term in eq. (|C,13j) . equalling — <7 2 C7p(l + 35/2)/4tt 2 
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T = 250 MeV 



T = 350 MeV 



T = 450 MeV 




Figure 8: The phenomenologically assembled scalar channel spectral function ps(u>), in units of 
w 2 Af 2 , for T = 250, 350, 450 MeV (from left to right). To the order considered, M is the heavy quark 
pole mass. Note that for better visibility, the axis ranges are different in the leftmost figure. As 
discussed after eq. (|C.14[) . we are not confident that these plots have a definite physical significance; 
the figures are meant for illustration only. 



times the leading term in eq. (|5.2|) . only by a multiplicative correction of the scalar density, 



^QCD — ^NRQCD 



fC F ( , 3 



1 - 2 — ^- 1 + -5 
8vr 2 V 2 



(C.14) 



Even though closer to unity than in eq. (|5.4|) . the normalization factor could be numerically 
significant. In fact, if we leave the normalization factor open, and search for a value mini- 
mizing the squared difference of the resummed and QCD results (with 5 = 0) in the range 
(oj — 2M)/M = 0.0 — 0.4 (thereby also accounting for thermal corrections), we find a best 
fit with an overall normalization factor 0.4 — 0.6, i.e. with a larger reduction than in the vec- 
tor case, in contrast to what eq. (|C.14p would sugg est@ This is perhaps another indication 
that the treatment of ps(u) within a Schrodinger-equation based resummed framework as in 
ref. [8] may not capture the correct physics. 

Nevertheless, putting this worry aside for a moment, we again construct an "assembled" 

i , (assembled) / (QCD) (resummed) \ mi -11 C j_1 1 • 

result as p s = max{ps ,p s ). ihe numerical value of the gauge coupling is 

taken from eq. (|5.5p . The outcome is shown in figs. [8] for (5 = and for various masses and 
temperatures, as a function of uj. Compared with the results in ref. [8], the overall magnitude 
is smaller by about 40 — 60%. At the same time, as is obvious from the plots, the two results 
do not interpolate to each other well; we have no explanation for this at the moment, but wish 
to repeat our concerns on the validity of the resummed near-threshold function ^ rcsummod )_ 



5 We note, however, that if we introduce another fit parameter, a horizontal energy shift, then the two 
results can be matched smoothly, with a multiplicative factor close to unity. We have not used this method 
in the plots because we prefer a universal procedure for the scalar and vector cases. 
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